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We investigate the possibility to identify photon-gluon fusion (PGF) events in polarised 
deep inelastic ep scattering, assuming the kinematics of the HERA collider, by a pair of 
■ charged high pt particles. In a Monte Carlo study we find possible selection criteria and 

CN . show the expected measurable asymmetries. We discuss the sensitivity to lS.G{x) and 

^ \ compare the result to the one obtained using di-jets to tag PGF events. 

1 Introduction 

On 

p • The extraction of the polarized gluon density AG is one of the future key measurements 
p^. in polarized physics. For the polarized HERA collider several methods have been proposed to 
D ', measure AG [jl|, the most promising one being the measurement of the di-jet asymmetry in deep 
^ ' inelastic scattering events 0]. In this paper we present a new method to access the polarized 
gluon density at HERA, using two charged particles with high transverse momentum pt, instead 
^ ■ of jets. This analysis follows closely the one proposed for the COMPASS experiment at 
CERN, detailed in ref. [Q. The technique to extract the (unpolarized) gluon density from high 
Pt hadrons has already been successfully applied at HERA, albeit based on one single high 
momentum charged track and for photoproduction data 0. 

The method is very similar to the one using di-jet events to extract AG. The aim is to isolate 
subprocesses which are sensitive to the gluon at the Born level. Fig. ^ shows that this is the case 
for the Leading Order (LO) photon gluon fusion process (PGF). In case of high partonic pt the 
two quarks will emerge as observable jets. These jets can be detected either calorimetrically, as 
exploited in p| or can be tagged via particles with high transverse momentum with respect to the 
7*1? axis in the 7*p centre of mass system (CMS), as studied here. The present experiments at 
the HERA collider are equipped with excellent tracking detectors in the photon fragmentation 
and central region, where most of the jets are expected. Hence we will use charged particles 
(tracks) for the high pt hadrons to tag the PGF process. 

In this exploratory study for HERA we are guided by the analysis of the di-jet events, and 
impose similar criteria to select the PGF process. Since in LO the two partons are produced 
with opposite pt, two charged hadrons with approximately opposite azimuthal angles in the 7*p 
CMS will be required. The background in LO results mainly from the QCD-Compton process 
(QCDC), see Fig. 1. For events with no hard QCD matrix element, namely 'quark parton 
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Figure 1: Feynman diagrams for the di-jet cross sections at LO, the Photon-Gluon Fusion (PGF) 
process (left) and the QCD-Compton process (QCDC) (right). 



model (QPM)' events, the pt of the hadrons only comes from fragmentation and intrinsic kt of 
the partons in the proton. Hence this background is expected to be strongly suppressed if two 
particles with a pt above 1-2 GeV are required. 

We make a full Monte Carlo simulation of the signal and background processes, include 
hadronization and higher order effects via parton showers. Starting from three different sets 
of polarised gluon distributions, shown in Fig. 2, we check the sensitivity of the measurements 
and extract AG{x). These distributions are the Gehrmann- Stirling (GS) sets A and C 0, which 
result from a QCD analysis of gi data, and the instanton-gluon distribution |^ . The latter results 
from a calculation of the polarised parton distribution in the Instanton Liquid Model [§. The 
distributions differ substantially and are purposely selected to demonstrate how poorly AG{x) 
is constrained by the present polarised data. For the unpolarised parton density functions the 
parametrizations of Gliick, Reya and Vogt in LO were used 0. 



2 Analysis method and expected asymmetries 



To estimate the sensitivity of high-pj track events to AG a study very similar to the one using 
di-jet events [0 was performed. The Monte Carlo program PEPSI 6.5 [|10], |Tl| was used and the 
kinematic range in Q'^ and y was taken to be: 5 < < 100 GeV^ and 0.3 < y < 0.85. PEPSI 
includes, apart from the polarized cross sections at LO, hadronization and (unpolarized) parton 
showers for higher order effects. The z — s scheme was used to regulate the divergencies of the 
matrix elements [0, |T^. 



For this analysis the acceptance of the HI detector [|1J] is assumed. Tracks are detected in 
either the forward or central tracking chambers. The acceptance for tracks in pseudorapidity 
7] in the laboratory frame amounts to the range: — 1.5 < rj < 2.5. Events were selected if two 
tracks were found inside this ?7-range each with a. pt > 1.5 GeV. The difference in the azimuthal 
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Figure 2: a) The true Sij distribution of the selected events (full line) and the Sij reconstructed 
from the two tracks (dashed line), b) The reconstructed Xg of high pt track events versus the 
true Xg of the event. For the reconstructed Xg a p^-dependent correction has been applied (see 
text). 



angle between the two tracks should be within 180 ± 60° and the difference in pseudorapidity 
should be \ritracki — Vtrack2\ < 2. In case the two highest pt tracks of the event do not fulfil these 
criteria, it was checked whether a third track with pt > 1.5 GeV was present, which does fit the 
criteria. In the following the convention will be that for the two selected tracks with Pi(l) and 
Pt(2) we have: Pt(l) > Pt(2). With these cuts and the parton distributions as given above the 
selected sample contains about 85% PGF events and 15% QCDC events. The QPM background 
is less than a percent. Events samples with a higher pt cutoff were studied as well, but the 
gain in purity and kinematic variable reconstruction (see below) was outweighted by the loss in 
statistics. 

A difference compared to the di-jet analysis is the reconstruction of the kinematics of the 
events. In the di-jet case an attempt is made to reconstruct the kinematics of the original parton. 
Thus, the invariant mass squared of the two jets is close to and well correlated with the true 
invariant mass of the hard subprocess. Here, where only two tracks opposite in azimuth are 
considered, the invariant mass squared s^J'^ = (Pi + ^2)^, with -Pi(2) being the four-momentum 
of track 1(2), is always much smaller than the true value of stj. This is shown in Fig. |]a. There 
is an offset between true and reconstructed Sij, but they are correlated. The correlation for 
Xg = x{l + Sij/Q"^) reconstructed from the tracks and the true value is shown in Fig. Here x 
is the Bjorken-x variable, and Xg is the momentum fraction of the proton carried by the gluon 
for the PGF process. The offset observed in Fig. ^ due to the incomplete Sij reconstruction is 
already corrected for in Fig. ^b. This was done by comparing the ratio of true Xg and Xg'^'^ in bins 
of Pt(2). A polynominal was fitted to logi^i^x^g'^'^) / log;^o(^9) ~ '^sP? (2) +aipt{2) +ao and the right 
side of this expression was used as a multiplicative correction factor to x^*^'^. The coefficients 
used were: = 0.855, ai = —0.052, 02 = 0.0018. Values of x^^^'^ presented in this paper are 
always corrected with this factor. It can be seen that the correlation is already rather good with 
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Figure 3: a) Measured asymmetries for high-pf track events, as a function of the true Xg. Asym- 
metries on the parton level (P.L.) are compared with the hadron level (H.L.) for different selection 
cuts, b) The measured asymmetry on hadron level versus the reconstructed Xg. For comparison 
the same parton level points of a) are shown again. The assumed luminosity is 200 pb~^. 



this simple correction factor. A more sophisticated analysis, e.g. using a real unfolding method, 
where bin-to-bin correlations are fully taken into account, can be performed when real data have 
been collected, and even better results can be expected. 

With the cuts detailed and used so far, the sample contains still a considerable amount, about 
23%, of events with Sij < 100 GeV^ for which NLO corrections could be potentially large (based 
on the experience with the di-jets). An additional criterium is used to reduce this fraction. One 
possibility is to increase the p^-requirement for the tracks. Raising it from 1.5 GeV to 2 GeV 
decreases the fraction of events with < 100 GeV^ to about 10%, but the total number of 
events is decreased by more than 40%. A better possibility seems to be to check the total 
transverse energy T^Et deposited in the calorimeter (i.e. with \ri\ < 2.8 in the laboratory frame) 
for the event. Requiring T,Et > 15 GeV reduces the fraction of low-Sjj events to about 10%, 
but removes only about 25% of the total statistics. For T,Et > 20 GeV a reduction of low-Sjj 
events to 5% can be obtained, while keeping 55% of the original statistics. The latter number is 
comparable to the result with the cut of pt > 2 GeV. We decide here for a cut on T,Et > 15 GeV 
as a reasonable compromise between statistics and low-Sjj contamination. 

The final event sample contains about 80,000 events/100 pb~^. The measurable asymmetries 
Ameas for the final event sample are shown in Fig. |^. Ameas is defined as in 0: 

Ameas i\^T I | _/\^TT ^^PpD A2tracks (1) 

The quantities N"^^ (A^TT) 

are the total number of observed events with two high pt tracks 
(A^Ti = Np\jp + N'q(jj-,q) with proton and electron spin antiparallel (parallel) to each other. 
A2tracks IS the truc physical asymmetry and on parton level the same as Adi-jet of ref. 0. The 
depolarisation factor D is given hj D = {y{2 — y))/{y'^ + 2(1 — y){l + R)), where R is the ratio 
of longitudinal to transverse 7*p cross section. 
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Figure ^ shows the asymmetries plotted versus the true Xg of the event. The error bars 
correspond to an integrated luminosity of 200 pb^^. The assumed polarised gluon distribution 
is gluon set A of Gehrmann and Stirling 0. Compared are the asymmetries on hadron level 
with and without the cut of T^Et > 15 GeV. The asymmetries are not changing significantly and 
the statistical errors for the low Xg bins are shghtly increased after the cut has been applied. 
Both asymmetries can also be compared to the asymmetry expected on parton level, where only 
kinematical selection cuts were applied plus a cut on the true Sij > 50 GeV^. The asymmetries 
on parton and hadron level agree very well, the statistical errors however increase when applying 
the selection on hadron level. In order to check that the asymmetry that can be related to AG 
we show in Fig. the asymmetry plotted as a function of the reconstructed Xg. The parton 
level asymmetry is exactly the same as in a), i.e. plotted versus the true Xg. Hence, also for 
reconstructed quantities an asjTiimetry is observed, not too different in size from the one at the 
parton level. 

The QCDC background is less than 15% in all bins except the two highest Xg bins where, with 
increasing Xg, the background increases to 30% and 50% respectively. The Aq/q distributions 
are however expected to be known with a precision of better than 10% at the time of this 
measurement, and can be subtracted. The background in the high Xg region could be also 
further suppressed selecting hadrons with opposite charge or with strangeness, as demonstrated 
ini- 



3 Extraction of AG 

The study of the extraction of AG/G and xAG follows very closely the procedure described 
in [0. We simulate 350 pb~^ of events with GS-A as input gluon density. This would in a real 
measurement correspond to the Monte Carlo generation of events and will therefore be called the 
'MC-set' here. Assuming that for each x-hin AG/G and the background corrected asymmetry 
Acorr are related to each other by a simple factor Ff. 




where i indicates the x-bin, we compute these factors using the MC-set. {A^orr = j^n+jvt'f ^ •) 
The factors were then multiplied with the asymmetry Acorr which corresponds to the measured 
asymmetry on hadron level in Fig. Three 'data sets' are considered: the standard one based 
on gluon set A, and two additional samples with gluon set C and the instanton induced 
gluon density [0 respectively. All of these data sets correspond to a luminosity of 200 pb~^. 
They were used to compute the asymmetries which were then multiplied with the same factors 
Fi calculated with gluon set A. The results are shown in Figs. ||a - It can be seen here that 
within the statistical accuracy the input distributions (solid lines) AG/G are properly extracted 
in all cases. A discrimination between different polarised gluon sets is thus possible. Figures |^d 
- ^ show the same result presented for the theoretically more interesting quantity xAG. The 
error bars are scaled errors of the left column, i.e. no uncertainty was assigned to G which is 
expected to be known to better than 5% at the time of this measurement |T^. 

The polarised gluon distributions used for the three data samples (see Fig. ^) differ sub- 
stantially but are all compatible with the present polarized data, stressing the need for direct 
measurements of AG{x). The GS-A and GS-C distribution show a similar small x behaviour, 
but differ considerably in the region around x ~ 0.1. The GS-C distribution is negative for this 
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Figure 4: For a luminosity of 200 pb ^ the sensitivity to extract AG/G (a-c) and xAG (d-f) is 
shown for different polarised gluon densities (see text). 



X region. The instanton-gluon is quite different from the GS sets. It remains negative over the 
full X range. The latter gluon is used in combination with the GS-A quark distributions for the 
study in this paper. 

The sensitivity to the shape of AG/G and xAG for an integrated luminosity of 500 pb^^ is 
shown in Fig. |^. The statistical errors for the x points are shown on the curves for AG/G (a-c) 
and xAG (d-f), showing the separation power of this measurement. 

The statistics is similar to the case when di-jet event are selected 0, in fact overall the 
track sample contains about 15% more events than the di-jet sample. However the events are 
differently distributed in the x^-bins. In this analysis we obtain more events at low Xg than in 
the di-jet case but less events at very high Xg. This is reflected in the statistical errors. The 
access to lower Sij values allows to explore the extraction of the polarised gluon distribution at 
lower Xg values, compared to the di-jet method. Due to the presently as yet unknown higher 
order corrections for the small Sij region, which could be potentially large, we have not pursued 
this further in this study. The rather big error in the highest Xg bin can probably reduced with 
a more sophisticated reconstruction method. With the simple correction method used in this 
analysis most events originating from this bin are smeared into the second high Xg bin. 

Comparing the di-jet and the 2-track sample event-by-event we find that about 40% of 
the 2-track events are also selected in the di-jet analysis. Thus a rather large fraction of the 
statistics is uncorrelated, i.e. we add statistically new and independent information to the di-jet 
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results. In addition the systematic errors for this measurement are partially different. This 
method is not sensitive to the calibration of the hadronic calorimeter and does not depend 
on a jet definition. One of the potentially most important systematics for this measurement 
is connected with the fragmentation functions. Therefore the analysis was repeated using the 
somewhat extreme independent fragmentation scheme instead of the LUND string fragmentation 
scheme in PEPSI. The results were found to be in good agreement with each other, showing 
that the result is not too sensitive to details of these functions. 



4 Conclusion 

The use of two high pt particles to extract the polarised gluon density at a polarised HERA 
has been studied. It has been found to have a similar potential to measure AG in the range 
0.002 < Xg < 0.2 as the di-jet measurement, with the additional possibility to reach even lower 
Xg values, if more sophisticated unfolding techniques can be used. At high Xg the precision of 
the di-jet method remains superior, though. Both methods are subject to different systematics 
and can therefore be considered as complementary. 
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